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Abstract 

We prove that the number of sphaleron instabilities of the 
ro-th SU(2) Einstein- Yang-Mills-Dilaton soliton and black 
hole equals n. 



1 Introduction 

The discovery of a countable family of solitons pQ and 
black holes j2] in static, spherically symmetric SU(2) 
Einstein- Yang-Mills (EYM) theory has triggered an ongo- 
ing interest in the study of matter models of this kind. The 
answer to the central question about the (linear) stability 
of these solutions turned out to be negative though 
The spherically symmetric perturbations which decouple 
into an even and an odd parity sector both contain nega- 
tive modes 0] E] ■ Whereas the former are called "grav- 
itational" (not having flat space-time analogues), the lat- 
ter allow for a "sphaleron-like" interpretation, cf. 0] |SJ E] 
and references therein. Numerical studies in |3j showed 
that the n-th soliton or black hole, characterized by the 
number n of nodes of the gauge potential, has exactly n 
negative gravitational modes. In the sphaleron sector, this 
has been shown numerically in 4 whereas in JSJEI the fol- 
lowing strategy has been used for an analytical proof: The 
arising system of coupled differential equations for the per- 
turbation amplitudes is cast into the form of a Schrodinger 
equation. Its formal supersymmetric partner Hamiltonian 
has a zero energy state which depends on an auxiliary 
function whose properties, due to the residual gauge free- 
dom of the perturbation equations, can be chosen such 
that the roots of the zero mode are the roots of the back- 
ground gauge potential. Since the supersymmetry trans- 



formation is reversible for strictly negative energies, the 
number of negative modes thus equals n. 

It is important to ask what role these EYM solitons and 
black holes play in string theory. One of the most widely 
studied stringy corrections stems from the superstring ef- 
fective action at low energies which supplements the EYM 
theory, among other scalar fields, by the massless neutral 
dilaton with exponential coupling to matter (Kaluza-Klein 
theories and cosmological models also suggest such a sup- 
plementation). This stringy generalization, the so-called 
Einstein- Yang-Mills-Dilaton (EYMD) theory, is specified 
by the action 
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where R is the scalar curvature, G Newton's constant, 
the dilaton field with coupling constant k, and F the 
Yang-Mills field with coupling constant A. After a suitable 
rescaling (cf. after ©), the theory is parametrized by 



7 = 



G 



(2) 
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The EYM and the Yang-Mills-Dilaton theory, both em- 
bedded in the EYMD theory, are recovered from the latter 
in the limits 7 — > and 7 — > 00, respectively. As for these 
extreme cases u\ , due to the breaking of scale invariance 
of the Yang-Mills action, a discrete familiy of static, spher- 
ically symmetric soliton and black hole solutions has been 
found in |S] for the intermediate values < 7 < 00 (for 
the extreme cases, rigorous existence proofs and classifi- 
cations of the solutions are given in 9 ). In this paper, we 
study the instabilities of these generalized solutions in the 
sphaleron sector. Using the strategy described above, we 
prove that the number of instabilities of the n-the soliton 
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or black hole equals the number n of nodes of the gauge 
potential. 

2 The number of instabilities 

In order to compute the action S from for the spheri- 
cally symmetric EYMD theory and, eventually, to derive 
the corresponding background and pulsation equations, we 
parametrize the spherically symmetric gravitational field 
by Schwarzschild coordinates, 



ds 2 = NS 2 dt 2 



dP 



r 2 {d9 2 + sin 2 9d(j) 2 ), 



(3) 



and the su(2)-valued gauge potential by 
A = aQT r dt + a\T r dr + (w — 1)(t^ d6 — tq sin 9 d<fi) 



-w(Tg d9 + T0 sin 9 d<fi), 



(4) 



cf. 10 . The metric coefficients N, S and the gauge 
amplitudes ao> <2i> w> w depend on r and t only. The 



spherical generators Tj of su(2), j 



are defined by 



Tj = r-ej, where r = a/(2i) with the usual Pauli matrices 
a = (<n, (72, (T3), and e r , eg = dge r , and = d^e r j sin# 
denote the unit vectors in the r, 0, and directions, re- 
spectively. 

Plugging the parametrizations © , into GJ > we compute 
the scalar curvature R using the Cartan structure equa- 
tions and the gauge field strength by F = dA+ [A,A]/2. 
Then, the action becomes 



If ( 1 e 2KV 
S = ~ drdt f -SX+S2 + —S3 



(•5) 



where the densities are given by the expressions 

Si = -S(N + rN'- 1), (6) 



S2 = r 2 I ^ - iVS^' 2 , 
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+ ATS ^ + a °*^ + ~ aow) 2 ) 
-ATS ((w' + aiw) 2 + (w' - aiw) 2 ) 

-^(- 2 +- 2 -l) 2 - W 

For the length S = VG/\, the rescaling of the form r t— * Sr, 
t 1— > t/<5, S 1— » <5 2 S, </? h- > ip/(\5), ciq 1— > <5ao, ai 1— * ai/5 
(iV, iw, w unchanged) allows to set k = 7 from J2J an d 
G = A = 1 to which choice we will stick in the following. 



In the static case, and with the magnetic ansatz ao = 
a\ = w = 0, the stationarity of the action (JSJ with ®;0> 
iJSJ , leads to the background Euler-Lagrange equations for 
the fields N, S, w, and ip, respectively, 



S' 


= s(rip' 2 +*e 2 ™w ,2 Y 


(9) 
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= 1 - r 2 Nip' 2 - U 


(10) 
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where we define 



U = 2e 2w [Nw' 2 + 



( w 2 - l) 2 
~2r~ 2 



(13) 



Since we are interested in a linear stability analysis for 
small spherically symmetric perturbations of the back- 
ground soliton and black hole solutions of i(§)l- i|13|) . we 
determine the pulsation equations for the perturbing fields 
5N, SS, Sao, 5a\, Sw, Sw, and Sip by linearizing the gen- 
eral Euler-Lagrange equations around these background 
solutions, cf. JT] for a general procedure. As for the 
EYM case mentioned in the introduction, the gravitational 
modes SN, 5S, Sw, Sip and the sphaleron modes Sag, Sa\, 
Sw decouple in the EYMD system in the magnetic ansatz 
do = a\ = w = 0. With the temporal gauge Sao = 0, the 
pulsation equations for the sphaleron modes which we are 
interested in become 



Sai = 
Sw = 

u Sw = 



2NS 

— - — (w Sw' — w'Sw — w 2 Sai) , (14) 
e -2iv NS ( e ^<P NS (Sw' - viSai))' 

^f - (w 2 - 1) Sw, (15) 



N 2 S 2 w'S ai 



(16) 



To simplify the form of these equations, we introduce the 
tortoise coordinate p defined by 



dp 
dr 
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(17) 



with p(0) = for the solitons, and p = —00 at the horizon 
of the black holes (from here on to the end of this section, 
the prime will denote the derivative with respect to p). 
Plugging the perturbations 



Sen =a{p)e lu \ 5w = f3{p)e i 



(18) 
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into the linearized dynamical equations (|14|) , (|15|) , , we 
are left with a coupled system of ordinary linear differen- 
tial equations for the gauge amplitudes a and (3 consisting 
of 
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5 / to an NS2 
-^{w'f3-w/3')-\ 5-u; 



(19) 
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+e 2 ~" ip NSw'a + (e 2j<p NSwa)' , (20) 
and the Gauss constraint 



erfwp = 0. 



(21) 



To rewrite the equations l|19|l , 1)20(1 , l|21|l in a more compact 
form, we define modified gauge amplitudes x, y, 
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y = e""p 



and modified background fields /, g, 
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Now, note that, if we multiply the right hand side of equa- 
tion l|19fl by e 2lv r 2 /S, take the derivative of the resulting 
product, and make use of the background equation (|11|) in 
the coordinate p, 
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(24) 



we get the right hand side of equation (|20|l multiplied by 
w. Hence, using the definitions l|22() . (|23|l . we can rewrite 



solve l]25p. In the following, we keep the Gauss constraint 
<|26[) also for the case u> — where it plays the role of a 
gauge fixing condition (also called the strong Gauss con- 
straint): The remaining gauge invariance in l(27|) is broken 
in the sense that h has to fulfill the linear second order 
differential equation 



9 . 



(28) 



Now, as in jSJEj, to determine the number of instabilities, 
we cast the system (I25II . H26(I into the form of a formal 
radial Schrodinger equation for the function (j> — x/f by 
eliminating the amplitude y with the help of the strong 
Gauss constraint We then get 



H°d> = uj 2 



(29) 



where the potential V° in the Hamiltonian H° = 
—d 2 /dp 2 + V° is given by 



= 2 



/'V ?'.f2 
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(30) 



Since V° has singularities at the roots of the gauge ampli- 
tude w = e~ lv f ', we try to construct the supersymmetric 
partner Hamiltonian H 1 of H° whose potential V 1 is such 
that the number of bound states of H 1 can easily be deter- 
mined. We proceed as it has been described in the intro- 
duction. Due to the residual gauge invariance (|2"7|) . . 
the Hamiltonian H° has a zero mode (j>o = h'/(fg). This 
zero mode allows for the facorization H° = Q + Q~ where 
the charges Q ± = ^fd/dp + a are expressed through the 
superpotential a = —cj)' /(f>o which solves the Riccati equa- 
tion V Q = g 2 — a' . Moreover, the potential of the super- 



equation JSJ as f 2 gx + f'y - fy' = iv 2 x, equation symmetric partner H 1 = Q~Q + = -d 2 /dp 2 + V 1 is given 



as (f 2 gx + f'y — fy')' — u> 2 fy, and the Gauss contraint 
(I21|l as lo(x' — fy) = 0. If u> ^ 0, the second equation is 
a consequence of the first and the third one. In this case, 
we are left with the system 



fy'-f'y = (f 2 g 
x' = fy- 



(25) 
(26) 



In the case w = 0, the Gauss constraint is void and we 
have to deal with the single pulsation equation ll'ol) which 
is invariant under the gauge transformation x > x + hi I g 
and y i— ► y + hf, where h is an arbitrary function. Hence, 



by V 1 = a 2 + a J . Using the definition z = —gh/h' and the 
strong Gauss constraint (|28|l , the superpotential takes the 
form cr = f / f + f 2 z. Therefore, since z solves the Riccati 
equation z' — f 2 z 2 — g, the first term on the right hand 
side of H30|) cancels. Hence, V 1 has the form 

V 1 = ^ + f 2 g + 2(f 2 z)'. 

Moreover, since /"// = NS 2 (w 2 - l)/r 2 + j(ip" + jtp' 2 ) 
from l|24|l. we get rid of the remaining singularities in V 1 
generated by w, 



ti 



y = hf 



(27) 



NS 2 



(3w 2 ~l)+ 1 ( V >" + Jv' 2 ) + 2(f 2 z)'. (31) 
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To find a candidate ipo for a zero mode of H 1 , we look for 
a solution of the equation Q + ipo = 0. Using the form of 
the superpotential a = f'/f + f 2 z from above, we have 



ipo = Cf exp dp f z 



(32) 



Let's now first treat the case of the solitons. Due to the 
remaining freedom in the choice of a solution of l|28[l . we 
construct a function z in Appendix IXlsuch that both V 1 is 
sufficiently regular and tpo is square integrable on the half- 
line and vanishes linearly at the origin: Plugging Ij38|) and 
(GHU into (J2II) and using the asymptotics of the background 
from Appendix [51 we find that V 1 is a bounded S'-wave 
potential with 



V 1 = C + 0(p), p^O, 
V 1 = —7 - - c j — — ] . p — > oo. 



P \ P 6 



(33) 



Furthermore, plugging Ij38(l . (|39|l into (|32|1 . the zero mode 
"00 behaves like 

O = Cp + 0(p 2 ), p^O, 

1 



0o = O 1-^1 , p ->• oo, 

which implies the desired properties of ipo (the non- 
uniqueness of V 1 does not affect the bound state prob- 
lem, cf. jHl). Now, due to Sturm-Liouville theory ^2Jj the 
number of bound states of H 1 equals the number of nodes 
of ipo- But, using definition l|23l) . the number of nodes of 
ipo equals the number of nodes of w because ipo is of the 
form H32|) . Since, as long as oj 2 < 0, the supersymmetry 
transformation is invertible, the claim follows. 

The derivation of the number of instabilities for the 
black hole case proceeds along the same lines as the one 
just presented for the soliton case, cf. Appendix IU1 for a 
brief discussion. 



3 Concluding remarks 

We studied the linear stability of the SU(2) EYMD soli- 
ton and black hole solutions of il§)l- H13l) under small time 
harmonic, spherically symmetric perturbations l|18[l with 
the help of the following strategy. The pulsation equa- 
tions (|14JI - (|16|I reduce to the Schrodinger eigenvalue prob- 
lem (I29|) with Schrodinger operator H° whose potential V° 
(|30|l has singularities located at the roots of the gauge am- 
plitude w. Due to the residual gauge invariance l|27|l . (|28|l . 



the Hamiltonian H° has a zero mode (po which allows for 
the factorization of H° in a product of the supersymmet- 
ric charges Q ± , cf. after H30|) . Moreover, using the gauge 
freedom in (|28J) . we can choose the potential V 1 H31(l of 
the supersymmetric partner Hamiltonian H 1 to be suf- 
ficiently regular. Since the zero mode ipo of H 1 has its 
zeroes at the nodes of w (KS2[| and since the supersymme- 
try transformation is invertible on the negative spectrum, 
the Sturm-Liouville theory implies that the number of in- 
stabilities equals the number of nodes of w. 

In J3| , for the EYM theory with a positive cosmological 
constant A, this strategy has been used to derive the linear 
instability of the sphaleron modes for all purely magnetic, 
static, spherically symmetric solutions, classified numeri- 
cally in Jl] and analytically in 15 with respect to the size 
of A and the number n of nodes of the gauge amplitude w 
(cf. as well |16j for solutions with similar properties of spa- 
tially compact sections in the EYM-Higgs theory) . In JJj , 
a numerical analysis in the gravitational sector revealed an 
unexpected dependence of the number of instabilities on 
the value of A: The number of unstable modes of the n = 3 
solution jumps from n to n — 2 as A crosses a critical value 
from below, a phenomenon which remains to be studied 
in more detail. 
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A Constructing z 

Using the definitions l(2*2*|l , (|2*3^l and the background equa- 
tion 10, we can rewrite (|28(l as a second order linear dif- 
ferential equation of Fuchsian type in the variable r, 

r 2 h" + 2r (l-w' 2 + ^-tp /2 +-frtp'^j ti - ^-h = 0. (34) 

This differential equation has two regular singular points, 
a first one at the origin, and a second one at infinity. As 
in we assume all the background fields to be analytic 
at the origin and at infinity (in the EYM case this has 
been proven in |S]). Hence, using the structure of the 
coefficients in Q34JI and the asymptotics from Appendix 151 
the Fuchsian theory yields the indices 1 and —2 at the 
origin, and 2 and —1 at infinity. Moreover, it turns out 
that the solution at the origin neither contains a 1/r nor 
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a logr term in its power series expansion. If we choose a 
special solution ho of 1(28(1 which vanishes at infinity, we 
have to keep the 1/r 2 term at the origin though, cf. [Sj. 
bincc, with 0, we assume continuity of the coefficients in 
iJsSJ, we are thus given a classical solution of 1(34(1 on the 
half-line with asymptotics 



h = £+0(1), 

C ( 1 
ho = - + -j 



0. 



(35) 
(36) 



To pick a function z having the desired properties de- 
scribed after 1(32(1 . we make use of the general solution 
of 1)28(1 constructed from the particular solution ho. In 
the coordinate p, the general z takes the form 



u 9 9 

' ho K~W 




(37) 



Plugging the expressions for the asymptotic behavior from 
Appendix 151 and ((35(1 . 1(36(1 into ((37(1 . we can achieve that 
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0{p), 



0. 



z = — - 



o 



p 



logp 



(38) 
(39) 



B Asymptotics of the background 

From |U we have the following asymptotic behavior of the 
fields at the origin (note that S can be expressed through 
w and (p with the help of the background equation Q), 

<p = (A) + 27e 2wo &V + 0(r 4 ), 

w = 1 - br 2 + 0(r 4 ), 

N = 1 -4e 2wo 6 2 r 2 + C(r 4 ), 

S = l + 4e 2wo 6 2 r 2 + 0(r 4 ). 

At infinity, the asymptotics looks like 
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2V 2 



Mr 



where M is the ADM mass and D stands for the dilaton 
charge. For a more detailed description and the meaning 
of the other constants, cf. 

C The black hole case 

The starting point is again equation 1(281) . In the black 
hole case, the coordinate p from ((17(1 runs over the whole 
real line with p = —oo at the horizon, r% say. The two 
regular singular points are now at p = ±oo. Rewriting 
equation ((34(1 in the coordinate q = (r — r/J/r^, assuming 
analyticity of the background solutions in some neighbor- 
hood of the horizon (cf. [5] for the EYM case), and 
using the asymptotics jS] 



w 
N 
S 



- +7t 
r 1 



V h 



V h 

w h (w 2 h - 1) 



q + 0(q 2 ) 
q + O (q 2 ) 



1-V h 
(1 - V h )q + O {q 2 ) , 

/ (2w 2 +j 2 V h )V h 



0(q 2 



where Vh = e 2lVh (w 2 — l) 2 /r^, the Fuchsian theory yields 
the indices 1 and (note that 1— Vh > 0). The asymptotics 
at infinity is described in Appendix[51with the indices from 
Appendix lAl Choosing the particular solution ho of (|2%f 
to be the power series belonging to index 1, we find, in the 
same way as in Appendix ^ that 



h 
h 



Cq 
Cq 



0(q 2 ), 
0(1), 



q^Q, 



(I 



Hence, with (|3*7|) , the general z has the form (|3*9l for p 
oo, and 



W 2 p 



P - 



Furthermore, the potential V 1 decays as 1/p 2 in leading 
order both for p — > oo and p — > — oo (cf. 133(1 ). and the 
zero mode ipo from 132(1 is square integrable on the whole 
real line, 



^o = 
^o = 



p^ -c 
p —>■ oo. 
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